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Precision studies of QCD at e + e~ colliders are based on measurements of event shapes 
and jet rates. To match the high experimental accuracy, theoretical predictions to 
next-to-next-to-leading order (NNLO) in QCD are needed for a reliable interpretation 
of the data. We report the first calculation of NNLO corrections (O(a^)) to three-jet 
production and related event shapes, and discuss their phenomenological impact. 

1 Introduction 

' Measurements at LEP and at earlier e + e~ colliders have helped to establish QCD as the 

theory of strong interactions by directly observing gluon radiation through three-jet pro- 
duction events. The LEP measurements of three-jet production and related event shape 
observables are of a very high statistical precision. The extraction of a s from these data 
sets relies on a comparison of the data with theoretical predictions. Comparing the different 
sources of error in this extraction, one finds that the experimental error is negligible com- 
pared to the theoretical uncertainty. There are two sources of theoretical uncertainty: the 
theoretical description of the parton-to-hadron transition (hadronisation uncertainty) and 
the uncertainty stemming from the truncation of the perturbative series at a certain order, 
as estimated by scale variations (perturbative or scale uncertainty). Although the precise 
■ size of the hadronisation uncertainty is debatable and perhaps often underestimated, it is 

certainly appropriate to consider the scale uncertainty as the dominant source of theoreti- 
cal error on the precise determination of a s from three-jet observables. From the planned 
luminosity of the ILC, one would expect measurements of event shapes comparable in sta- 
tistical quality to what was obtained at LEP, thus allowing for precision QCD studies at 
ILC energies. 

So far the three-jet rate and related event shapes have been calculated [1,2] up to the 
next-to-leading order (NLO), improved by a resummation of leading and sublcading infrared 
logarithms [3,4] and by the inclusion of power corrections [5]. 

QCD studies of event shape observables at LEP [6] are based around the use of NLO 
parton- level event generator programs [7]. As expected, the current error on a s from these 
observables [8] is dominated by the theoretical uncertainty. Clearly, to improve the de- 
termination of a s , the calculation of the NNLO corrections to these observables becomes 
mandatory. We present here the first NNLO calculation of three-jet production and related 
event shape variables. 
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2 Calculation 



Three-jet production at tree-level is induced by the decay of a virtual photon (or other 
neutral gauge boson) into a quark- antiquark-gluon final state. At higher orders, this process 
receives corrections from extra real or virtual particles. The individual partonic channels 
that contribute through to NNLO are shown in Table [1] All of the tree-level and loop 
amplitudes associated with these channels are known in the literature [9-12]. 

For a given partonic final state, jets are reconstructed according to the same definition 
as in the experiment, which is applied to partons instead of hadrons. At leading order, 
all three final state partons must be well separated from each other. At NLO, up to four 
partons can be present in the final state, two of which can be clustered together, whereas at 
NNLO, the final state can consist of up to five partons, such that as many as three partons 
can be clustered together. The more partons in the final state, the better one expects the 
matching between theory and experiment to be. 

The two-loop 7* — > qqg matrix elements 
were derived in [9] by reducing all relevant 
Feynman integrals to a small set of master 
integrals using integration-by-parts [13] and 
Lorentz invariance [14] identities, solved 
with the Laporta algorithm [15]. The mas- 
ter integrals [16] were computed from their 
differential equations [14] and expressed 
analytically in terms of one- and two- 
dimensional harmonic polylogarithms [17]. 

The one-loop four-parton matrix ele- 
ments relevant here [11] were originally de- 
rived in the context of NLO corrections 

m i i -, n . , ., ,. ,, . , to four-jet production and related event 

lablc 1: Partonic contributions to thrce-ict , ' ,„, „ . , . 

„ , . . shapes 18,19 . One of these four-iet parton- 

hnal states m pcrturbativc OOD. , , , , . , , 

level event generator programs [19] is the 

starting point for our calculation, since it 
already contains all relevant four-parton and five-parton matrix elements. 

The four-parton and five-parton contributions to three-jet-like final states at NNLO con- 
tain infrared real radiation singularities, which have to be extracted and combined with the 
infrared singularities [20] present in the virtual thrcc-parton and four-parton contributions 
to yield a finite result. In our case, this is accomplished by introducing subtraction func- 
tions, which account for the infrared real radiation singularities, and are sufficiently simple 
to be integrated analytically. Schematically, this subtraction reads: 
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where dcr NNLO denotes the real radiation subtraction term coinciding with the five-parton 
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tree level cross section da^NLO in all singular limits [21]. Likewise, d<r]^£ is the one- loop 
virtual subtraction term coinciding with the one-loop four-parton cross section der]^ io in 
all singular limits [22] . Finally, the two-loop correction to the three-parton cross section is 
denoted by dcr]^ L0 . With these, each line in the above equation is individually infrared 
finite, and can be integrated numerically. 

Systematic methods to derive and integrate subtraction terms were available in the liter- 
ature only to NLO [23,24], with extension to NNLO in special cases [25]. In the context of 
this project, we fully developed an NNLO subtraction formalism [26], based on the antenna 
subtraction method originally proposed at NLO [19,24]. The basic idea of the antenna 
subtraction approach is to construct the subtraction terms from antenna functions. Each 
antenna function encapsulates all singular limits due to the emission of one or two unresolved 
partons between two colour-connected hard partons. This construction exploits the univer- 
sal factorisation of phase space and squared matrix elements in all unresolved limits. The 
individual antenna functions are obtained by normalising three-parton and four-parton tree- 
level matrix elements and three-parton one-loop matrix elements to the corresponding two- 
parton tree-level matrix elements. Three different types of antenna functions are required, 
corresponding to the different pairs of hard partons forming the antenna: quark-antiquark, 
quark-gluon and gluon-gluon antenna functions. All these can be derived systematically 
from matrix elements [27] for physical processes. 

The factorisation of the final state phase space into antenna phase space and hard phase 
space requires a mapping of the antenna momenta onto reduced hard momenta. We use 
the mapping derived in [28] for the three-parton and four-parton antenna functions. To 
extract the infrared poles of the subtraction terms, the antenna functions must be integrated 
analytically over the appropriate antenna phase spaces, which is done by reduction [29] to 
known phase space master integrals [30]. 

A detailed description of the calculation will be given elsewhere [31]. 

3 Results 

The resulting numerical programme, EERAD3, yields the full kinematical information on a 
given multi-parton final state. It can thus be used to compute any infrared-safe observable 
in e + e~ annihilation related to three-particle final states at 0(a 3 s ). As a first application, 
we derived results for the NNLO corrections to the thrust distribution [32] . 

In the numerical evaluation, we use M z = 91.1876 GeV and a s (M z ) = 0.1189 [8]. 
Figure [T] displays the perturbative expression for the thrust distribution at LO, NLO and 
NNLO, evaluated for LEP and ILC energies. The error band indicates the variation of the 
prediction under shifts of the renormalisation scale in the range /i £ [Q/2;2Q] around the 
e + e~ centre-of-mass energy Q. 

It can be seen that even at linear collider energies, inclusion of the NNLO corrections 
enhances the thrust distribution by around 10% over the range 0.03 < (1 — T) < 0.33, 
where relative scale uncertainty is reduced by about 30% between NLO and NNLO. Outside 
this range, one does not expect the perturbative fixed-order prediction to yield reliable 
results. For (1 — T) — > 0, the convergence of the perturbative series is spoilt by powers of 
logarithms ln(l — T) appearing in higher perturbative orders, thus necessitating an all-order 
resummation of these logarithmic terms [3,4], and a matching of fixed-order and resummed 
predictions [33]. 



LCWS/ILC2007 




The perturbative parton-level prediction is compared with the hadron-level data from 
the ALEPH collaboration [34] in Figure [TJ Similar data are also available from the other 
LEP experiments [35] . The shape and normalisation of the parton-level NNLO prediction 
agrees better with the data than at NLO. We also see that the NNLO corrections account 
for approximately half of the difference between the parton-level NLO prediction and the 
hadron-level data. A full study including resummation of infrared logarithms and hadroni- 
sation corrections is underway. 

4 Conclusions 

We developed a numerical programme which can compute any infrared-safe observable 
through to which we applied here to determine the NNLO corrections to the thrust 

distribution. These corrections are moderate, indicating the convergence of the perturbative 
expansion. Their inclusion results in a considerable reduction of the theoretical error on the 
thrust distribution. Our results will allow a significantly improved determination of the 
strong coupling constant from jet observables from existing LEP data as well as from future 
ILC data. 
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